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ABSTRACT 

Context. The metal content of dwarf galaxies and the metal enrichment of the intergalactic medium both suggest that mass loss 
from galaxies is a significant factor for the chemical evolution history of galaxies, in particular of dwarf galaxies. However, no clear 
evidence of a blow-away in local dwarf galaxies has been found so far. 

Aims. Dwarf galaxies provide a perfect environment for studying feedback processes because their kinematics and their generally 
low gravitational potential support the long-term survival of shells, filaments, and holes. We therefore performed a detailed kinematic 
analysis of the neutral and ionised gas in the nearby star-forming irregular dwarf galaxy NGC 2366 in order to make predictions about 
the fate of the gas and to get a more complete picture of this galaxy. 

Methods. A deep Ho- image and Fabry-Perot interferometric data of NGC 2366 were obtained. They were complemented by H i 
synthesis data from the THINGS survey. We searched for line-splitting both in Ha and H i by performing a Gaussian decomposition. To 
get an idea whether the expansion velocities are high enough for a gas blow-away, we used the pseudo-isothermal halo model, which 
gives us realistic values for the escape velocities of NGC 2366. The good data quality also allowed us to discuss some peculiarities of 
the morphology and the dynamics in NGC 2366. 

Results. A large red-shifted outflow north west of the giant extragalactic Hn region with an expansion velocity of up to 50kms~' is 
found in Ha, but not in Hi. Additionally, a blue-shifted component north of the giant extragalactic Hn region was detected both in 
Ha and Hi with an expansion velocity of up to 30kms~'. A comparison with the escape velocities of NGC 2366 reveals that the gas 
does not have enough kinetic energy to leave the gravitational potential. 

Conclusions. This result is in good agreement with hydrodynamic simulations and suggests that we need to examine even less massive 
galaxies (M m = 10 6 M Q ). 
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1. Introduction 

Irregular dwarf galaxies are known to be the sites of giant star 
formation regions. The feedback between massive stars and the 
interstellar medium (ISM), and on larger scales also the inter- 
galactic medium (IGM), is one of the most important processes 
in the evolution of these galaxies. Massive stars are strong 
sources of radiation and mechanical energy. Photoionisation is 
the most probable mechanism, but also shocks that are produced 
by stellar winds and supernova (SN) explosions inject huge 
amounts of energy into the ISM. This leads to numerous ionised 
gas structures in and around the galactic plane, visible on 
deep Ha images (e.g.,|Hunter & Gallagher] 1997||Martin| 1998 



Bomans et al. 1997; van Ey meren et al.|2007] >. However, ionised 
gas also exists at kpc distances away from any place of current 
star formation ( Hunter et al.| 1993 1. These structures might 
be fragmented cool (T * 10 4 K) shell structures that were left 



behind by the expanding hot gas (Mac Low & Ferrara 1999), 
but they can also be explained by turbulent mixing between the 
hot and the neutral gas (Slavi n et al.|1993[ ). 



In order to explain the observations, theoretical models were 
developed in which the gas is shock-heated by collective SNe 
and accelerated into the ISM forming a thin shell of swept-up 
ambient gas. Because of Rayleigh-Taylor instabilities the shell 
can rupture and the hot gas is expelled through tunnel-like 
features, called chimneys, into the halo of the host galaxy 



( Norman & Ikeuchi|1989| l. On the way into the halo and in the 
halo itself, the gas cools down and depending on the strength of 
the gravitational potential, the hot gas might fall back onto the 
galactic disc in the shape of cool clouds, which is described in 
the Galactic Fountain scenario ( Sh apiro & Field|1976[ ). 
In some cases, the gas can be accelerated to velocities beyond 
the escape velocity of the host galaxy. Especially the relatively 
low escape velocities of dwarf galaxies will facilitate the 
removal of substantial amounts of interstellar matter. Therefore, 
the question comes up whether the gas stays gravitationally 
bound to the galactic disc (outflow) or whether it can escape 
from the gravitational potential by becoming a freely flowing 



wind (galactic wind). Numerical simulations by |Mac Low & 
|Ferrara| (jT999j model superbubble blowout and blow-away in 
dark matter dominated dwarf galaxies of different mass and with 
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Fig. 1. /?-band image (left panel) and continuum-subtracted Ha image (right panel) of NGC 2366. The optical centre is indicated by 
a white cross on the /?-band image. The contrast of the Ha image is chosen in a way to emphasise the small-scale structure. The 
supergiant shell SGS 1 (see |van Eymeren et al.|2007| l as well as two of the main H n regions are marked. In order to stress weaker 
structures and to differentiate them from the noise, we used adaptive filters based on the H-transform algorithm (Richter et al. 1991). 



different energy input. Their models show that only in low mass 
dwarfs (M gas = 10 6 M ) a significant fraction of the hot gas can 
escape from the gravitational potential. The metals however, 
produced by massive stars and released during a SN explosion, 
have a high probability to be blown away. The fraction is almost 
independent of the gas mass or the energy input of the host 
galaxy. 

Therefore, star-forming irregular dwarf galaxies provide a 
perfect environment to study feedback processes and to hunt for 
galactic winds. A number of studies already concentrated on the 



lar (IB(s)m) dwarf galaxy ( |de V aucouleu rset al.| 1991 . We 
adopt a distance of 3.44 Mpc from Tolstoy et al. ( 1995 i that 



kinematics of ionised gas structures (e.g., Hunter & Gallagher 
1997 1 |Martin|[l99"8"} |van Eymeren et al.||2007) , o thers on the 



distribution and kinematics of the neutral gas (e.g., |Thuan et al.| 
2004| >. Only a few studies address both components together in a 
detailed analysis. Additionally, most authors use single long-slit 
spectra for the examination of the ionised gas component, which 
usually only cover some parts of extended target objects. During 
the last decade, a lot of progress has been made in the field of 
the so-called 3d spectroscopy. This makes it possible to observe 
large parts of a galaxy or even the whole galaxy together with 



sufficient spectral information in one exposure (e.g., Wilcots & 
|Thurow|2001|). 



In this paper we present a study of the neutral and ionised gas 
component in the nearby irregular dwarf galaxy NGC 2366. 
The Ha emission as a tracer of the ionised gas component was 
observed with a Fabry-Perot (FP) interferometer that provides us 
with a complete spatial coverage of the stellar disc and relevant 
spectral information. The FP Ha data are complemented by 
optical imaging and H i synthesis data. 

NGC 2366 is classified as a barred Magellanic-type irregu- 



places NGC 2366 in the M8 1 group. The nearest neighbour is 
NGC 2403 at a projected distance of 290 kpc, which makes any 
kind of recent interaction unlikely. Its appearance in Ha is domi- 
nated by the Giant Extragalactic H n Region (GEHR) NGC 2363 
in the south-west at 07 h 28 m 29.6 s , +69 d ll m 34 s (see Fig.[T] right 
panel) with a luminosity twice as bright a s 30Doradus (|Chu 
|& Kennicutt|19 94). Recent observations by van Eymeren et al. 
( |2007| l revealed the existence of numerous ionised gas structures 
up to kpc-size, especially close to the GEHR where most of the 
current star formation activity takes place. Several of them have 
been shown to expand. The expansion velocities, however, were 
in all cases much lower than the escape velocities of the galaxy. 

This paper is organised as follows: The observations and the 
data reduction are described in § 2. Sect. 3 is a comparative de- 
scription of the Ha and H i morphology. § 4 presents a kinematic 
analysis of both Ha and H i data. The results are subsequently 
discussed in § 5, followed by a short summary in Sect. 6. 

2. Observations and data reduction 

2.1. Optical imaging 

A 900 s /?-band image and a deep - 3600 s - Ha image of 
NGC 2366 were obtained on the 16th of October 2007 with the 
Calar Alto 3.5m telescope equipped with MOSCA, the Multi 
Object Spectrograph for Calar Alto. After the standard data re- 
duction performed with the software package IRAF, the contin- 
uum image was subtracted to produce an image of the pure Ha 
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line emission. In order to emphasise weaker structures and to 
differentiate them from the noise, we used adaptive filters based 
on the H-transform algorithm ( Richter et al. 1 1 99 1\ . Both images 
are shown in Fig.[T] 

2.2. The Fabry-Perot Ha data 

FP interferometry of NGC 2366 was performed on the 1st of 
March 2006 with the 1.93m telescope at the Observatoire de 
Haute-Provence, France. We used the Marseille's scanning FP 
and the new photon counting camera ( |Gach et al.| |2002). The 
field of view is 5?8x5'8 on the 512x512 pixels of the detec- 
tor and is slightly limited by the interference filter to 5'.5 x5'5, 
which gives a spatial resolution of O'.^S per pixel. The Ha 
line was observed through an interference filter centred at the 
galaxy's rest wavelength of 6564.53 A with a FWHM of 10 A. 
The free spectral range of the interferometer - 376 km s -1 - was 
scanned through 24 channels with a sampling step of 15 kms -1 . 
The final spectral resolution as measured from the night sky lines 
is about 50 kms -1 . The seeing was between 3" and 4". 
60 cycles were observed with an integration time of 10 sec per 
channel, hence 240 sec per cycle. After removing bad cycles, a 
total integration time of 232 min remained. We used a neon lamp 
for the phase and the wavelength calibration. The data reduc- 
tion was done with the software package ADHOCwQ written by 
Jacques Boulesteix. 

2.3. The Hi data 

We were provided with a fully-reduced data cube from "The 
Hi Nearby Galaxy Survey" (THINGS, [Walter et al.||2008l >, a 
high spectral and spatial resolution survey of Hi emission in 
34 nearby galaxies obtained with the VLA in its B, C, and D 
configuration. The synthesised beam size is 13" x 12" after us- 
ing a "natural" weighting. The spectral resolution is 2.6 kms -1 . 
For more details see Walter et al. (2008 1. We applied a 3-point 
Hanning smoothing to the cube, which improved the noise level 
from 0.6 to 0.4 mJy beam -1 . 

The moment maps were created from the un-smoothed cube by 
removing all emission below a 2.5 cr threshold. The processing 
and the subsequent analysis of the H i data was performed with 
GIPS^ 



3. General morphology 

3. 1 . Results from optical imaging 

In addition to the /?-band and the Ha images presented in Fig. [T] 
a second grey-scale presentation of the Ha image with a differ- 
ent contrast in order to emphasise the weaker structures is given 
in Fig. 



A.l 



The galaxy has a bar-like appearance. The main star- 
forming complex is located at the southern end, coinciding with 
the GEHR NGC 2363 (Fig. [T] right panel). The H n region west 
of the GEHR is sometimes referred to as a satellite galaxy (e.g., 
Driss en et al.|20 00). Whether an independent system or not, this 
H ii region is connected to the stellar disc by diffuse, filamentary 
gas structures. 

A catalogue of ionised gas structures can be found in 



Eymeren et al. (2007 1. As the field of view of the Ha image 
used by van Eymeren et al. (2007} only covered the main part 



of the stellar disc, we made some new detections on the larger 



image used here. To the very north of the disc a very interesting 
structure was found, which is the superbubble SB4 (continuing 
the labelling of van Eymeren et al.|2007[ ). At the same position, 
the /?-band image shows a small cluster of stars. This bubble is 
remarkable because its morphology consists of a sharply edged 
shell pointing to the north and a smooth shell pointing south. 
Additionally, there is a second shell, SB5, heading south-east 
at a distance of 1 kpc from SB4 without an ionising star cluster 
nearby. Another peculiarity in NGC 2366 are the very small and 
faint H n regions that are scattered along the disc on both eastern 



and western side with no apparent connection to it (see Fig. A. 1 



3.2. Results from Hi synthesis data 

Figures [2] and [3] show the H i channel maps of the cube af- 
ter applying a 3-point Hanning smoothing, superimposed on a 
greyscale presentation of our /?-band image. The white cross in 
the first channel marks the optical centre of the galaxy. The beam 
is placed into the lower left corner of the same channel. The cor- 
responding heliocentric velocities are printed in the upper right 
comer of each channel. The Hi emission is very extended (up to 
a factor of two in comparison to the /?-band image) and diffuse. 
The most prominent structure is a tail-like feature at a velocity 
of 100 km s -1 . 

The H i moment maps as well as the global intensity profile are 
displayed in Fig. |4] The Hi intensity distribution (upper left 
panel, 0th moment map) looks very patchy with several inten- 
sity maxima in the inner parts and a more diffuse filamentary 
structure in the outer parts. On a larger scale, the Hi forms two 
bright elongated structures or ridges that are embedded in fainter 



and smoother gas (see also Hunter et al. 2001 Thuan et al. 
2004). Both ridges are parallel to the major axis and the west- 
ern one coincides with the tail-like feature detected in the chan- 
nel maps. The total Hi flux density is 213 Jy km s -1 , which is in 



1 http://www.oamp.fr/adhoc/ 

2 The Groningen Image Processing System 



good agreement with single dish measurements by, e.g., |Thuan 
|& Martin| ( |1981| l. This shows that the D array of the VLA is very 
sensitive to extended emission. The Hi mass derived from the 
0th moment map is given in Table [T] 

3.3. A comparison of the neutral and ionised gas distribution 

Figure [5] shows a comparison of the ionised and neutral gas dis- 
tribution in NGC 2366. In the upper panel, the Hi contours are 
overlaid over the continuum-subtracted Ha image. The optical 
galaxy lies in the centre of the Hi distribution, i.e., in the area 
of highest H i intensity. The H i maximum flux density coincides 
with the GEHR, the region of currently the strongest star forma- 
tion in NGC 2366. 

For a more detailed study, Fig. [5] lower panel shows an enlarge- 
ment of the main optical part of NGC 2366. In order to better 
differentiate between H i maxima and minima, the H i contour at 
0. 1 Jy beam -1 is overlaid in black, the H i intensity contours from 
0.15 to 0.4 Jy beam -1 in steps of 0.05 Jy beam -1 are overlaid in 
white. First, it can be seen that the H i maximum coinciding with 
the GEHR is slightly offset from the centre of the Ha emission. 
The shift is of the order of one beam size or even more and there- 
fore not a resolution effect. Such offsets have also been found in 
other galaxies like, e.g., IC 10 ( |Hodge et al.||1994| ) and are ex- 
plained by sequential star formation. 

Additionally, the arm-like feature in the north-east as well as the 
western H n region coincide with an H i intensity maximum. And 
finally, several smaller Hi intensity peaks can be found in the 
northernmost part of NGC 2366, close to but offset from SB4. 
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Fig. 2. H i channel maps of NGC 2366 (contours) as obtained from the VLA using "natural" weighting of these data, superposed 
on the /?-band image. A 3-point Hanning smoothing was applied to the cube improving the noise level from 0.6 mJy beam" 1 
to 0.4 mJy beam" 1 . The original channel spacing is 2.6 km s -1 . Contours are drawn at -1.2 (-3<x), 1.2 (3<x), 3, 6, 12, and 
24mJybeam" 1 . The synthesised beam is displayed in the lower left corner of the first channel map. The optical centre of the 
galaxy is marked by a white cross in the same channel map. The corresponding heliocentric velocities are printed in the upper right 
corner of each channel. 
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SB5 seems to be located in an Hi hole, which gives rise to the 
assumption that it expands into the neutral medium by work- 
ing like a snowplough. The neutral gas in front of the shell is 
compressed and maybe shock-heated. In the back, a cavity of 
low densities evolves visible as an Hi hole. This phenomenon 
cannot be observed in the vicinity of SB4, but the filamentary 
structures north west of the GEHR are all found in an area of 
low Hi column density (see Sect. 



5.1 



4. Kinematic analysis 

In this section, we will analyse the kinematics of NGC 2366 in 
order to study the characteristics of the gas and to explain the 
peculiarities in its morphology. Spectra extracted from the Ha 
and the H i data cubes show that both emission lines are some- 
times split into several components. Tasks which create veloc- 
ity fields from a 3d data cube cannot take line-splitting into ac- 
count. Instead, line-splitting can then even lead to an incorrect 
measurement of the line properties. Therefore, we performed a 
Gaussian decomposition by interactively fitting the Ha and H i 
emission (IRAF task splot) in order to measure the characteris- 
tics of the lines, especially the peak velocity, and to create the 
velocity fields from these values. Only detections above a 3<x 
limit were considered. All given velocities are heliocentric ve- 
locities measured along the line of sight. 

4. 1 . The Ha velocity field 

The Ha velocity field (Fig. |6j upper panel) shows the compo- 
nent of the strongest intensity. In order to match the seeing and 
to improve the signal to noise ratio, we summed over 3x3 pixels 
before performing the Gaussian decomposition. The overall ve- 
locity gradient runs from the south-west with velocities of about 
65kms~ 1 to the north-east with velocities of about 120 km s -1 . 
We found two major expanding gas structures (marked in black), 



a large red-shifted and a faint blue-shifted outflow. 
On the lower panel of Fig. [6] we present as an example one of 
the spectra extracted from the FP data cube (black solid line). 
The Ha line is split into two components, one at 70kms~ I (blue 
(dark grey) long-dashed line) and one at 130 km s _1 (red (light 
grey) long-dashed line). The sum of both Gaussian fits is plotted 
with a green (light grey) short-dashed line and is in good agree- 
ment with the observed spectrum. This is one of the spectra ex- 
tracted from the area of the large red-shifted outflow (marked by 
a cross in Fig. [6] upper panel). As can be seen, the outflow domi- 
nates the Ha profile, whereas the main component is only visible 
as a wing on the left side. The lines to the left side of the Ha line 
are night sky lines which could not properly be removed. 
This outflow has been detected befor e by|Roy et al. ( 1991)) also 
using a FP, and by Martin] ( fl 998} and [van Eymeren et al. \ 2007) 
performing long-slit echelle spectroscopy. |Roy et al.|( 1991 1 con- 
cen trated their analysis on the GEHR NGC 2363, |Martin| ( |1998) l 
and van Eymeren et al.] ( 2007|) were limited by the slit size. E.g., 
only one slit position by van Eymeren et ak] ([2007) intersects 
the outflow, which made it impossible to define its whole ex- 
tend. This new FP exposure now shows the complete size of the 
outflow. Assuming a distance of 3.44 Mpc, we estimated a to- 
tal length of 1 .4 kpc, which is two times greater than measured 
by |van Eymeren et ak ([2007 ) and even four times greater than 
observed by Roy et al. ( 1991 ). This length makes it one of the 
largest outflows ever detected in a dwarf galaxy. The gas expands 
with a velocity of up to 50kms _I red-shifted in comparison to 
the rotational gradient, which is almost a factor of 2 higher than 
the detections of Roy et ak]( |1991|l and van Eymeren et aL] ( |200"7) . 
In case of the observations by Roy et al. ( 1991 1, the difference in 
velocity is most probably due to the fact that they only detected 
the edge of the outflow where the velocities are lower due to its 



symmetry. The slit position of van Eymeren et aL] ( |2007[ ) inter- 
sects the whole outflow. Here, the discrepancy could be due to 
their very limited field of view, which did not give them a true 
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Fig. 4. The Hi moment maps of NGC 2366 using "natural" weighting, which leads to a synthesised beam of 13" xl2". Top left: The 
Hi intensity distribution (Oth moment). Contours are drawn at 0.005, 0.01, 0.025, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7 Jy beam" 1 
where 5mJybeam _1 correspond to a column density of 3.6 x 10 19 atoms crrT 2 . Top right: The Hi velocity field (1st moment). 
Contours are drawn from 40 to 160 km s _1 in steps of 10 km s _I . The systemic velocity of 97 km s _I is marked in bold. Bottom left: 
The velocity dispersion (2nd moment), overlaid are the same Hi intensity contours as on the Oth moment map. Bottom right: The 
global intensity profile of the galaxy. The short-dashed line marks zero intensity. 



Ha rotation curve, or due to their imprecise H i velocity mea- 
surements, which they used as a reference value for outflowing 
gas. The comparison with the H i data presented here will show 
that the expansion velocity is indeed higher than measured so far 
(see Sect.|4~4l>. 

The other detected outflow is a faint blue-shifted component 
north of the GEHR with velocities of 60 to 70kms _1 in com- 
parison to the general velocities of 90 to 100 kms -1 , which is 
discussed in more detail in Sect. 15.11 

As the Ha line is only split in two areas across the galaxy and 
as the outflow is both times the dominant component, only the 
main Ha velocity field is shown here. 



4.2. The Hi velocity field 

In order to analyse the general Hi kinematics, we worked on 
the moment maps created with GIPSY (see Sect. 12.31. The 



moment maps created with GIPSY (see Sect. 
Hi velocity field (1st moment map), displayed on the upper 
right panel of Fig. [4] is fairly regular in the inner parts, but 



disturbed in the outer parts. The velocity gradient goes from the 
south-west with velocities of about 60kms _1 to the north-east 
with velocities of about 145kms _1 , which is similar to the Ha 
velocity gradient. The isovelocity contours close at both ends of 
the galaxy, which might be caused by the bar and possibly by 
a warp structure in the outer part. Especially the north-western 
part shows a completely different kinematic behaviour to what 
is expected from a regularly rotating galaxy. 
The velocity dispersion map (Fig. [4j lower left panel, 2nd 
moment) varies between lOkms 1 in the outer parts and 
17kms _1 in the inner parts of NGC 2366. The dispersion 
peaks with values of up to 28kms~ 1 in an area east of the 
optical disc which is also clearly offset from any Hi maxi- 
mum. Close to the Ha shell SB4 in the north, a very small 
Hi dispersion peak can be seen with a maximum value of 
25kms~ 1 , which is again offset from a nearby Hi maximum. 
This is not unusual as a high velocity dispersion suggests a 
high level of turbulence which will naturally decrease the Hi 
column densities. The high Hi velocity dispersions on the 
western and eastern edges of the galaxy have to be treated 
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Fig. 5. A comparison of the Ho- and Hi morphology. The up- 
per panel shows the continuum-subtracted Ha image. Overlaid 
in black are the H i intensity contours at 0.005, 0.01, 0.025, 0.05, 
0.1, 0.2, 0.3, and 0.4 Jy beam 1 . The lower panel displays an en- 
largement of the optical part. The Hi contour at 0.1 Jy beam 1 
is overlaid in black, the Hi intensity contours from 0.15 to 
0.4JybeatrT 1 in steps of 0.05 JybearrT 1 are overlaid in white, 
which allows us to distinguish between maxima and minima in 
the neutral gas distribution. 



with care. A closer look at the Hi velocities in these areas will 



show that this is an effect of low signal to noise (see Section 4.4 1 



4.3. The Hi rotation curve 

In order to investigate the fate of the gas, we need to know some 
of the kinematic parameters of NGC 2366 like its inclination or 
its rotation velocity. Therefore, we derived a rotation curve from 
the H i data by fitting a tilted-ring model to the observed velocity 
field. Figure [7] upper left panel shows the resulting curves. 
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Fig. 6. The Ha velocity field of NGC 2366. The upper panel 
shows the velocity field of the strongest component. The posi- 
tions of the two detected outflows are marked. The lower panel 
displays one example spectrum (black solid line) which was ex- 
tracted from the FP data cube at the position of the huge red- 
shifted outflow (black cross in the Ho- velocity field) and which 
shows the flux in arbitrary units vs. the velocity offset (not cor- 
rected for the systemic velocity of NGC 2366). The Ha line 
is split into a component at 70kms~' (blue (dark grey) long- 
dashed line) and one at 130kms _I (red (light grey) long-dashed 
line). The sum of both Gaussian fits is plotted with a green (light 
grey) short-dashed line. 



At the beginning, initial estimates for the kinematic parameters 
had to be defined, which were obtained by interactively fitting 
ellipses to the Hi intensity distribution using the GIPSY task 
ellfit. These were then used as an input for the tilted-ring fitting 
routine rotcur. The width of the rings was chosen to be half the 
spatial resolution, i.e., 6" in this case. In order to get the most 
precise values, three different approaches were made by always 
combining receding and approaching side. First, the initial es- 
timates were all kept fixed. The resulting curve is indicated by 
the green (light grey) symbols in Fig. [7] upper left panel. In a 
second approach, the parameters were iteratively defined for all 
rings. Up to a radius of 300", which corresponds to a distance 
of 5 kpc from the dynamic centre, no significant deviation or a 
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Fig. 7. The Hi kinematics of NGC 2366. Top left: The rotation curve derived from performing a tilted-ring analysis. Different 
approaches were chosen: the black symbols represent the best-fitting parameters, the error bars indicate receding and approaching 
side, respectively. The green (light grey) curve was derived by taking the initial estimates and keep them fixed, the red (dark grey) 
curve by taking the best-fitting parameters and let them vary. Top right: The observed H i velocity field. Bottom left: The model 
velocity field, based on the best-fitting parameters. Bottom right: The residual map after subtracting the model from the original 
velocity map. The outermost Hi intensity contour at 0.005 JybearrT 1 is overlaid in black. 



sudden change of any of the parameters was noticed so that an 
average value (given in Table [TJ was taken for each parameter 
(black symbols). In this case, a rotation curve was also measured 
for receding and approaching side alone, indicated by the error 
bars of the black symbols. As a last approach, the so derived pa- 
rameters were all kept free in order to reproduce the result of the 
second approach (red (dark grey) symbols). 
The green (light grey) symbols are in very good agreement with 
the black ones over the whole range of 450 ". The red (dark grey) 
symbols follow the black and green (light grey) ones up to a ra- 
dius of 200". Beyond a radius of 300", the deviation is of the 
order of 30 to 40kms~' and the differences between receding 
and approaching side become larger. The reason for this is that 
the filling factors of the rings drop from about 1 to about 0.5 at 
a radius of 300" as the galaxy does not have a clear elliptical 
shape. This leads to a higher uncertainty in calculating the rota- 
tion velocity. Therefore, all values above this radius have to be 
treated with care. 

In general, the shape of the derived rotation curve fits the ones 
expected for a dwarf galaxy. In the inner 100 " the velocity gradi- 
ent is very steep and linear, indicating solid body rotation, which 
is a characteristic sign of dwarf galaxies. From 100 " to 200 " the 



curve is still linear, but less steep coming to a plateau at 220". 
From 320" on, the rotation curve declines, which has already 
been implied by the closing velocity contours in the 1st moment 
map (Fig. |4] upper right panel). 

The best-fittingnarameters derived in the iterative approach are 
given in TableTT] The systemic velocity of 97kms~ 1 is in good 
agreement w i th the values measured by Hunter et al.| ( 2001 1 and 



Thuan et al. 



2004| i (99 and 101 km s respectively). We de 
rived an inclination of 64° and a position angle of 45° in com- 
parison to i — 65° and PA = 46° published by, e.g., Hunter et al. 
(2001 1. Our derived rotation velocity of 50km s is also in very 
good a greement with the values measured, e.g., by|Hunter et al. 



(2001) of 47 km s" 1 . 
In order to prove the reliability of the derived parameters, a 
model velocity field with the best-fitting parameters was created 
using the GIPSY task velfi (Fig. [7] lower left panel) and sub- 
tracted from the original velocity map (upper right panel). The 
residual map can be seen in Fig. [7] lower right panel. Overlaid in 
black is the outermost H i intensity contour. The overall velocity 
field is well represented by our derived parameters except for the 
outer parts in the north-west and south-east. Here, the residuals 
reach values of more than 20 km s _I in comparison to a general 
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Table 1. Hi properties measured from the THINGS data. 



Parameters [Unit] 


NGC 2366 


optical centre": 




a (J2000.0) 


Q7 h 2gm 54 gs 


6 (J2000.0) 


+69° 12' 57" 


D [Mpc] b 


3.44 


dynamic centre c : 




a (J2000.0) 


07 h 28 m 55.4 s 


6 (J2000.0) 


+69° 12' 27" 


v sys [kms~'] c 


97 


■ r]c 


64 


PA [°T 


45 


v rot [kms-'] c 


50 


F m [Jy kms 1 ] 


213 


M ffl [10 s Mo] 


5.95 


H i diameter ['] 


15x7 


[kpc] 


15x7 


H i / opt. ratio 


1.9x2.8 


< cr > [kms -1 ] 


10/ 17 


crp e ak [kms~'] 


28.2 


>HI,max [kpc] 


7.5 


M dyn [10" M ] 


4.3 



a Data from NED. 

b Distance measured from cepheids l Tolstoy et al.|l995 



Derived by fitting a tilted-ring model to the H i data 
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Fig. 8. A comparison of the neutral and ionised gas kinematics. 
The residuals after subtracting the H i velocity field from the Ha 
velocity field are shown. 



He velocity field. Therefore, the FP data were smoothed to fit 
the Hi spatial resolution of 13" x 12". The residual map is 
shown in Fig. [8] At most positions, the velocities of the neutral 
and ionised gas are in good agreement with offsets of less than 
± 10 km s _1 . The large red-shifted outflow clearly stands out with 
an expansion velocity of up to 50 km s . Additionally, the whole 
northern part except for the arm-like feature to the east seems to 
have a faint blue-shifted component with an offset maximum of 
30kms _1 and a median offset of 15kms _1 . This expanding gas 
could already be seen on the original Ha velocity map. It coin- 
cides with an area of diffuse and filamentary emission between 
the GEHR and the northern tail (see Fig. [5] lower panel). For a 
detailed catalogue of the single filaments see |van Eymeren et al.| 
( |2007| , their Fig. A. 1 and Table Al . 

In order to look closer at the H i kinematics in the area of the ex- 
panding ionised gas and to check the regions of high H i velocity 
dispersion, we performed a Gaussian decomposition of the Hi 
data as described above. The result is shown in Fig. [9] and some 
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example spectra extracted from the H i cube are given in Fig 
together with the fitted Gaussian profiles for the single compo 
nents (long-dashed blue (dark grey) and red (light grey) lines) 
and the resulting sum (short-dashed green (light grey) lines). 
The velocities were averaged over one beam size, which is in 
this case roughly 12" in both spatial directions. The strongest 
component is presented in the middle panel. Overlaid in white 
are the Hi velocity dispersion contours at 20 and 25kms~ 1 and 
the outer Ha intensity contour in black. For a comparison, the 
blue- and red-shifted components are shown on the left and right 
panel. Note that regions where we did not find a blue- or red- 
shifted component were filled with the main component. 
Fig. 1 1 shows an enlargement of SB4 and the area of the GEHR. 
Overlaid over the three Hi velocity components are some of 
the Ha intensity contours. The superbubble SB4, the supergiant 
shell SGS1 (see Fig. [TJ, and the outflow OF are marked. As al- 
ready mentioned in Sect. |3.3| the huge outflow north west of the 
GEHR expands into an area of low column density (a factor of 
at least two lower than the surroundings). The Hi velocity maps 
in Fig. 1 1 lower row show no line-splitting, which can also be 
seen by looking at a spectrum in this area (Fig. 10 panel d). Only 
close to the origin of the outflow, a second component appears 
in the H i line profile. With a velocity of 100 km s _1 in compari- 



it is most probably not 



son to the outflow velocity of 130 km s _1 , 
connected to the outflow. 

SB4 has also no counterpart in H i (see Fig. 1 1 upper row). The 
same is true for SGS 1 (lower row). Line-splitting of the H i emis- 
sion takes place in the north-eastern part of the GEHR going to 
the south, which is at the position of the ridge mentioned by 
|Hunter et aLl ( |200T) or[Thuan et al.| ( |2004] >. This is discussed in 
more detail in Sect. 15. 31 

Finally, a hint of the faint blue-shifted component detected on 
the Ha velocity map (see Fig. [6]) is also visible in Hi with the 
same velocity of 65 kms -1 (see Fig. [To] panel c). 



value of ± lOkms The north-western part could not be fitted 
with this model. 

4.4. A comparison of the neutral and ionised gas kinematics 

The Ha velocity field (see Fig. [6]l shows two major deviations 
from the overall rotation velocity, a red-shifted component north 
west of the GEHR and a blue-shifted one north of the GEHR. 
As a next step, the Hi velocity field was subtracted from the 



5. Discussion 

Our analysis has shown that NGC 2366 harbours many interest- 
ing features both in its gas distribution and its gas kinematics. In 
the following subsections, we try to explain the detections (and 
also a non-detection) and deal with the important question what 
the fate of the outflows is. 
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Fig. 9. Gaussian decomposition of the H i line. Blue-shifted (left panel), main (middle panel) and red-shifted (right panel) compo- 
nents of the H i velocities are shown. Overlaid in white are the H i velocity dispersion contours at 20 and 25 km s _1 together with the 
outermost Ha intensity contour in black. The positions of the four spectra shown in the next figure are indicated. 
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Fig. 11. The same as Fig.^for two peculiar areas in the central region of NGC 2366. The upper row shows the superbubble SB4 and 
its surroundings, the lower row presents the GEHR and the onset of the northern tail. Again, the three Hi velocity components are 



plotted. Overlaid in black are the Ha intensity contours. The superbubble SB4, the supergiant shell SGS1 (see van Eymeren et al. 
|2"007] > and the outflow OF are marked. 
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Fig. 10. Some examples of Hi line profiles. The Gaussian pro- 
files fitted to the single components are plotted with blue (dark 
grey) long-dashed lines and red (light grey) long-dashed lines, 
the resulting sum of the profiles is displayed with a green (light 
grey) short-dashed line, (a) The eastern ridge with an additional 
red-shifted component, (b) The western ridge with an additional 
blue-shifted component, (c) A hint of the blue-shifted compo- 
nent north of the GEHR with Vheiio — 65kms~ 1 . (d) The red- 
shifted outflow north-west of the GEHR with Vheiio — 130km s _I 
is not detected in H i. 



5.1. The outflows 

Two major outflows were found in the Ha velocity field. The 
large red-shifted outflow expanding from the GEHR to the north- 
west is not detected on the H i velocity map, whereas the blue- 
shifted outflow north of the GEHR has a counterpart in Hi. 
Additionally, we see shell-like structures (SGS1, SB4) on the Ha 
image with no kinematic evidence on the Hi map. Resolution 
effects cannot be responsible for the non-detection as the struc- 
tures cover a field of several beam sizes. Unfortunately, we can- 
not trace the shells in the Ha velocity field as SB4 lies outside the 
field of view and SGS1 is covered by artificial emission caused 



in the detector (see Sect. 5.2 1. However, we have already shown 
that SGS1 is kinematically visible on long-slit echelle spectra 
( |van Eymeren et al.|2007| l. 

The question which now arises is why the blue-shifted outflow is 
detected both in Ha and H i, whereas the red-shifted outflow is 
only detected in Ha. This is probably an effect of age and energy 
input. Considering their position in the galaxy, it can be assumed 
that both outflows have a different origin. The large red-shifted 
outflow is most probably driven by the two large star clusters 
in the GEHR ( |Drissen et al.|200l) . Assuming a deprojected dis- 
tance from the two central star clusters of approximately 1 kpc 
and a constant expansion velocity of 50kms _1 , we can estimate 
the expansion time of the red-shifted outflow to be 2 x 10 7 yr. 
According to Driss en et al.| ( |200T) , the older star cluster has an 
age of 3 to 5 Myr. This means that the outflow must originate 
from a former star formation event and is now driven further out 
by the new event. The estimated cooling times for both outflows 
(see Martin| 1998[ ) lie clearly above the expansion ages so that ra- 
diative losses can be neglected. Therefore, the gas could already 
have been fully ionised by the former event. On the other hand, 
it is also possible that our assumption of a constant expansion 
is wrong. The expansion into a dense environment could have 
caused a decline in the expansion velocity of the gas over the ex- 
pansion time. A continued decline would of course have strong 
consequences for the fate of the gas. 

The blue-shifted gas probably gets its energy from a star cluster 
north of the GEHR. Assuming a distance of 500 pc and again 
a constant expansion velocity, this time of 30kms _1 , the ex- 
pansion time is comparable to the one of the huge red-shifted 
outflow. As the age of this star cluster is not known, we cannot 
tell whether it is alone responsible for the blue-shifted outflow. 
What we can tell is that because of the comparable expansion 
ages, the energy input for the blue-shifted component must have 
been lower or the neutral gas densities higher, as the gas is not 
fully ionised, yet. Furthermore, the fact that we can only detect 
a blue-shifted component does not necessarily exclude the ex- 
istence of a red-shifted component. A detection depends on the 
column density. In this case, it could be possible that the blue- 
shifted gas runs into an area of higher column-density, whereas 
the red-shifted part expands into an area of lower column den- 
sity. The same could be true for the large red-shifted outflow. 
Here, we do not see a blue-shifted component in Ha and we do 
not see any counterpart in H i. Possibly, the energy input was so 
strong that all the gas was immediately ionised, e.g., by shock 
heating. The expansion velocities are almost twice the velocities 
of the blue-shifted outflow. The low H i column density in com- 
parison to its surroundings also suggests that high amounts of 
neutral gas have already been ionised. 

Both processes, ionisation of the neutral gas and blowout should 
affect a smooth Hi intensity distribution. Figure 12 shows the 
central region of the 1st moment map. As already mentioned in 
Sect. |3.2| the intensity distribution is quite patchy. At several 
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interpretation as an expanding superbubble is indeed intimating. 
Therefore, we suppose that most probably, Roy et al. ( 1991 1 in- 
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Fig. 12. The locations of the central H i holes in NGC 2366. The 
inner part of the H i intensity distribution is shown. In order to 
facilitate the comparison of the locations of the H i holes (indi- 
cated by red (light grey) ellipses) with the ones of the outflows, 
the outer Ha contour is displayed in black. 



locations associated with filamentary and shell-like structure in 
Ha, we see small holes in the Hi distribution indicated by red 
(light grey) ellipses. These Hi holes are often close to the ob- 
served outflows. Especially the large red-shifted outflow is ex- 
panding into an area of very low H i column density. 
In order to complete the discussion of the gas kinematics and 
the morphology, the hot ionised gas has to be included, too. 
For NGC 2366 ROSAT and XMM Newton data are available. 



Stevens & Strickland| ([1998 ) analysed the ROSAT data which 
show that if there is X-ray emission associated with NGC 2366, 
it is barely detectable. The same is true for the XMM image 
where we might see faint X-ray emission in the centre of the 
GEHR (which is to be expected as it is the main area of current 
star formation activity) as well as in the northern part. However, 
no strong source associated with NGC 2366 could be detected. 

5.2. Superbubble blowout in NGC 2366? 



As already mentioned in Sect. 4.1 Roy et al. (1991 1 observed 
NGC 2366 with a FP interferometer centred on [Om] /I5007. 
Their velocity field shows line-splitting only within the GEHR 
NGC 2363. Next to the huge outflow they also reported the de- 
tection of an expanding superbubble around NGC 2363 with 
a diameter of 200 pc and an expansion velocity of 45kms~ 1 . 
Echelle spectroscopy of the central region performed by Martin 
( 1998| and |van Eymeren et al.| ( |2007] ) revealed neither a blue- nor 
a red-shifted component close to the centre of NGC 2363. The 



bubble is also not visible in the H i velocity field (see Fig. 1 1 



During the analysis of our FP data, we looked very carefully at 
the GEHR. What we found is that not only in this but also in 
other galaxies, the FP observations of bright emission line re- 
gions, e.g., of GEHRs, seem to be affected by blurring inside the 
detector. This leads to an unusually high, artificial velocity gra- 
dient in these emission line regions. In NGC 2366, e.g., the blur- 
ring caused a velocity gradient of more than 200km s within 
1 kpc. This is unphysical in comparison to the slowly rising ve- 
locity gradient in the other parts of the galaxy and is not observed 
in H i. Assuming that this gradient comes from true emission, its 



terpreted an artefact as a real expanding bubble as they were not 
aware of the technical problem of blurring inside the detector. 
They could also have been biased by a saturation problem of the 
photon counting system, leading to an artificial splitting of the 
observed emission line. 



5.3. Spiral arm structure in NGC 2366 

NGC 2366 shows several peculiarities. Photometric HST obser- 
vations by Tikhonov & Galazutdinova| ( 2008 1 reveal an overden- 
sity of blue stars at a galactocentric distance of 0.9 kpc against 
the overall decrease in the young-star number density with galac- 
tocentric distance, which has also been suggested by a B/V pre- 
sentation by |Hunter & Elmegreen ( 2006| l. At the same distance, 
our Hi observations as well as the ones by |Hunter et al.| (2001) 
and |Thuan et al.| ( [2004[ ) show the presence of two parallel ridges 
running along the major axis. Another peculiarity is the H n com- 
plex west of the GEHR, which is often referred to as a satellite 
galaxy interacting with NGC 2366 (e.g., |Drissen et al.|2000 1. 
Tikhonov & Galazutdinova ( 2008 ) suggest from their results of 
the stellar photometry that the H i ridges, interpreted as a depro- 
jected Hi ring by |Hunter et al.| ( |2001| i and |Thuan et al.| ( |2004) , 
could be weak spiral arms. Our multi-wavelength study gives 
further evidence for the existence of two weak spiral arms in 
NGC 2366, one located at the eastern part of the stellar disc and 
one at the western side (see Fig. |A. 1[ >. In Ha, the spiral arms can 
be traced by an alignment of small isolated H n regions at large 
distances from the disc. The most distinct ones (numbered 7 and 



8) are 2.7 and 3 kpc away from the disc (see also |Hunter et al. 
[2001) . 

Looking again at the enlarged view of the central region (Fig.[5]l, 
we see maxima in the H i column density that follow the spiral 
arms out of the disc. This could indeed be explained as a con- 
centration of neutral gas formed by a density wave going through 
the galaxy. All along the spiral arms, the H i line is split into two 
components, which has already been mentioned by, e.g.JThuan 



et al. (2004 1. A closer look at the Hi velocities (Fig.|9]l shows 



that next to the main component, a red-shifted component runs 
along the eastern arm to the south and a blue-shifted compo- 
nent along the western arm to the north (see also the example 
spectra in Fig. [TUj panels a and b). The eastern arm crosses the 
supergiant shell SGS 1 . As the red-shifted component is visible 
on large scales, it is unlikely that we can see an additional red- 
shifted component belonging to the SGS 1 . 
Furthermore, we compared our data to the GALEX FUV im- 
age of NGC 2366 ( |Gil de Paz & Madore||2005) where the star 
distribution is more pronounced than on the Ha image. It also 
suggests the existence of two weak spiral arms, which coincide 
with the H i ridges and the small H n regions. 
Taking all this together, we can confirm the conclusion of 
|Tikhonov & Galazutdinoval (|2008;) that NGC 2366 seems to be a 
weak two-armed spiral. 



5.4. Outflow or galactic wind? 

We found two main expanding gas structures in NGC 2366. One 
is only visible in Ha with an expansion velocity of 50kms _1 , 
the other one shows similar expansion velocities in Ha and H i 
of 30 km s~ 1 . Both were probably driven by a former star forma- 
tion event. 

We now want to make some statements about the fate of the gas 
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by comparing the expansion velocities of the outflows to the es- 
cape velocities of NGC 2366. As dwarf galaxies are dark matter 
dominated, we can calculate the escape velocities from dark mat- 
ter halo models. There are two main competing models which 
describe the density distribution (and therefore the appearance of 
the rotation curve) in a dark matter dominated galaxy, the NFW 
model derived from Cold Dark Matter simulations (e.g., [Navarro] 
|et al.||1996"t |Moore et al~||1998| |1999| l which predicts a cuspy 
core, and the empirically derived pseudo-isothermal (ISO) halo 
( |Binney & Trem aine 1987 ) predicting a core of constant density. 
This disagreement between simulations and observations, the so- 
called cusp-core discrepancy has been discussed in many publi- 
cations during the last decade (e.g., |de Blok & Bosma 2002). 
Recent studies have shown that the rotation curves of dwarf and 
low surface brightness galaxies can indeed better be described 
by an ISO halo, at least in the inner few kpcs (Kuzio de Naray 
|et al.|2008||Spano et al.|2008) . 

As the detected outflows are all closely located to the dynamical 
centre of NGC 2366, we decided to model the escape velocities 
by using the ISO halo. This decision was further supported by 
a mass decomposition which confirms that the rotation curve of 
NGC 2366 can better be described by the ISO halo ( |van Eymeren 
2008). Its density profile is given by 



PisoO) = po 



1 + 



2\ 



(1) 



with po being the central density and r c the core radius. The es- 
cape velocity is then given by 



= V 2 4 +iog (^)) 



(2) 



with v c being the circular velocity and r max being the maximum 



radius of the dark matter halo (see Binney & Tremaine|1987| l. 
We plotted the resulting escape velocities in Fig.|13| The Hi 
rotation curve including receding and approaching side is indi- 
cated by small grey triangular-shaped symbols. The expanding 
gas structures are marked by black large triangles. We corrected 
the values for an inclination of 64° as calculated in Sect. |4.2| 
which leads to an increase in velocity of about 10%. The circular 
velocity was measured from the rotation curve to be 50kms _1 
(see Table [TJ. We plotted the escape velocity for two different 
radii r max = 7.5 kpc (dotted line) and r max = 15kpc (solid line). 
The lower value of r max was chosen to equal the size of the H i 
outer radius. This is the outermost point where we can directly 
measure a contribution of dark matter from the rotation curve, 
which gives us therefore a lower limit for the radius of the dark 
matter halo. Most probably, the dark matter halo is much larger 
so that we chose the second value to be twice the size of the H i 
radius (which might still be too small). Nevertheless, a higher 
r max increases the escape velocities, which means that all radii 
greater than the H i radius decrease the probability of a galactic 
wind so that the Hi radius can indeed be regarded as a lower 
limit for the escape velocity. 

Fig. [13] shows that both outflows have expansion velocities that 
lie clearly below the escape velocity. As the ISO halo is charac- 
terised by a core of constant density, the escape velocities in the 
inner parts are even higher than by using a rapidly decreasing 
NFW density profile. This result is in good agreement with hy- 
drodynamic simulations (Mac Lo w & Ferrara|1999] see Sect.[TJ. 
With a gas mass of about W M Q NGC 2366 is l ocated at the 
upper end of the mass distribution simulated by Mac Low &| 
|Ferrara| (fl999). This means that according to their simulations, a 
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Fig. 13. Escape velocity for a pseudo-isothermal halo of r max = 
7.5 kpc (dotted line) and r max = 15 kpc (solid line). The ob- 
served rotation curve is indicated by small grey triangles. The 
error bars represent receding and approaching side. The expand- 
ing gas structures are marked by black large triangles. 



blow-away of gas is very unlikely. Nevertheless, the escape ve- 
locity drops significantly at greater distances from the dynamic 
centre, which enhances the probability of a galactic wind for gas 
in the halo. However, the deep Ha image of NGC 2366 shows no 
filamentary or shell-like structures at several kpcs from the disc. 
Either the emission is too weak to be detected or the gas has not 
managed so far to expand to such large distances. 

6. Summary 

Optical images, Fabry-Perot interferometric data, and Hi 
synthesis observations were used to get new insights into the 
morphology and the kinematics of the neutral and ionised gas 
components in the nearby irregular dwarf galaxy NGC 2366. 
The most important results are shortly summarised here. 

In agreement with recent studies we suggest the existence 
of two weak spiral arms in NGC 2366. 

The galaxy harbours two major outflows detected in the Ha 
velocity field. The large red-shifted outflow has an expansion 
velocity of 50kms~'. The FP data allowed us for the first time 
to measure the whole extent of this outflow. With a length of 
1 .4 kpc it is one of the largest outflows found so far in a dwarf 
galaxy. This outflow as well as the supergiant shell SGS1 and 
the superbubble SB4 do not have any counterpart in H i and they 
all expand into an area of low H i column density. We therefore 
suggest that major parts of the gas have already been ionised. 
The second outflow is located north of the GEHR NGC 2363 
and expands with 30kms _1 blue-shifted. It is also visible in Hi. 
Several nearby small H i holes show that a fraction of the neutral 
gas has already been ionised or blown away. An estimation 
of the expansion ages reveals that the current star clusters are 
probably not the first drivers of the outflows or that the gas is 
slowed down which decreases the chance for a galactic wind. 
We compared the measured expansion velocities with the escape 
velocities of the galaxy, calculated from a pseudo-isothermal 
halo model. In both cases and independent of the choice of 
r max , the expansion velocities of the outflows stay far below the 
escape velocity, which means that the gas is still gravitationally 
bound and would need an event like a SN explosion to be ac- 
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celerated to the escape velocity. Our result is in agreement with 
hydrodynamic simulations and draws therefore our attention to 
less massive galaxies. 



of the potential spiral arms are indicated. Additionally, an exten- 
sion of the catalogue by van Eymeren et al. ( 2007 1 can be found. 
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Table A.l. The most prominent structures and their sizes in 
NGC 2366 - an extension -. 



Source 


Diameter 


Diameter v hc i io FWHM 




["] 


[pc] [kms -1 ] [kms -1 ] 


SB4 


30 


500 


SB5 


25 


417 
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Appendix A: Ha image and extension of the 
catalogue of ionised gas structures 

Here, the continuum-subtracted Ha image is again presented in 
an enlarged version and with a different contrast as in Fig. [T] 
to emphasise the faint and diffuse ionised gas structure (see 
Fig. |A. 1 [ >. The isolated Hn regions as well as the new detections 
of shell-like material are marked and numbered. The positions 
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Fig. A.l. Continuum-subtracted Ha image of NGC 2366. The contrast is chosen in a way to demonstrate the large-scale structures. 
The extraplanar H n regions are marked by white boxes and numbered in black. Additionally, the two northern superbubbles are 
numbered following van Eymeren et al. (2007 1. The positions of the spiral arms are indicated by white lines. 



